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In a variety of rare-earth based compounds singlet superconductivity coexists with helical mag- 
netism. Here we demonstrate that surfaces of thèse Systems should generically host a finite density 
of zero-energy Majorana modes. In the limit of vanishing disorder, thèse modes lead to a divergent 
contribution to zero-energy density of states and to zero-temperature entropy proportional to the 
sample surface area. When confined to a wire geometry, a discrète number of Majorana modes can 
be isolated. The relative ly large characteristic energy scales for superconductivity and magnetism, 
compared to other proposais, as well as the lack of need for fine-tuning, make helical magnetic su- 
perconducting compounds favorable for the observation and expérimental investigation of Majorana 
fermions. 
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I. INTRODUCTION 

Majorana fermions (MF) have attracted considérable 
attention recently as a promising tool for topological 
quantum Computing, quantum information storage^^^ 
and for unconventional quantum transport phemonensP. 
The appeal of MF is in their non-local character, when 
realized in topological states of matter: a single binary 
quantum degree of freedom can be split into two MF 
spatially isolated from each other, and thus protected 
from decoherence caused by local environmental pertur- 
bations. An active search for material realizations of 
MF is under way. A list of the candidate Systems in- 
cludes chiral p-wave superconductors^ ^, topological in- 
sulator with proximity induced superconductivit}'^, non- 
centrosymmetric superconductoré^^, and 5/2 Quantum 
Hall stat^. In thèse Systems Majorana zéro modes can 
appear in the vortex cores or at the edge&mtl^. 

The practical implementation of most of thèse propos- 
ais is hampered by the small energy scale which sépa- 
râtes the Majorana modes from the electronic contin- 
uum, and by the need to precisely tune différent System 
parameters. For instance, recently, an élégant proposai 
to create and manipulât e Majorana fermions was made 
based on superconducting nano-wires in the présence of 
Rashba spin-orbit interaction and magnetic field^^ It 
requires a délicate balancing of the varions coupling con- 
stants: a gap in the electronic spectrum, which has to 
be sufEciently large to prevent detrimental effects of dis- 
order, has to be opened by an external magnetic field, 
which, at the same time, cannot be too large to avoid 
destroying superconductivity. It is clear that progress 
in the expérimental search of Majorana fermions would 
be greatly facilitated by the identification of candidate 
Systems, where Majorana physics can be explored in a 
larger energy range, without the need to satisfy highly 
demanding expérimental constraints. 

Here, we propose the use of helical magnetic supercon- 
ductors (HMS) as such candidate System. While the bulk 



properties of HMS have been investigated in the past, the 
sub-gap surface states in this class of materials have re- 
mained unexplored. We show that the spin-helix in thèse 
Systems plays a rôle analogous to the combination of the 
spin-orbit interaction and of the magnetic field in the case 
of superconducting Rashba nanowires, and leads to for- 
mation of zero-energy Majorana surface states. Impor- 
tantly, however, the pair-breaking effects of the exchange 
interaction in HMS are dramatically reduced due to the 
spatial modulation of magnetism, as is evidenced by a 
large variety of existing materials^^^ that exhibit coex- 
isting helical magnetism and superconductivity with rel- 
atively large transition températures (e.g, in ErNi2B2C 
the superconducting critical température is 10.5 K and 
the Neel température is 6.8 K). This property makes the 
readily available HMS a favorable Systems for the détec- 
tion of Majorana zéro modes, without the need to match 
demanding expérimental constraints, and offering an in- 
crease of an order of magnitude - or more - in the en- 
ergy scale relevant for Majorana physics as compared to 
a typical semiconductor based System. While in ideal- 
ized truly ID Systems the HMS and the superconducting 
Rashba Systems are unitarily équivalent, in wires of finite 
width or in bulk Systems - relevant to actual experiments 
- they are not. As we demonstrate, a remarkable feature 
of 2D and 3D HMS is that they generically possess a fi- 
nite density of zero-energy Majorana surface modes. The 
resulting large density of states can be accessed with sur- 
face probes (tunneling microscopy or photoemission) and 
should also manifest in thermodynamical measurements 
as a residual low-temperature entropy proportional to the 
surface area of the sample. 

The interplay of superconductivity and magnetism has 
long been a subject of active research. Uniform ferromag- 
netism is known to be strongly antagonistic to singlet su- 
perconductivity, due to the orbital magnetic field, which 
disrupts the superconducting phase by introducing vor- 
tices, and the exchange field that tends to split the singlet 
Cooper pairs. However, if the magnetization is spatially 
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non-uniform on a scale smaller than the superconducting 
cohérence length, both of thèse effects are dramaticany 
suppressed. In fact, the onset of superconductivity itself 
can drive a uniform ferromagnet into a non-uniform (he- 
hcal, or Ising-domain) state^^ Examples of materials 
where hehcal magnetism coexists with superconductivity 
are those where magnetism originates from the partiahy 
fihed f-orbitals of rare earth atoms, as is the case for the 
compounds such as HoMoeSg, ErRh4B4, TmNi2B2C^. 
Because the f-orbitals are very compact, they directly in- 
teract only with the itinérant électrons, which can dérive 
from the higher energy delocalized bands of the same 
atoms, or from other éléments in the compound. The 
exchange interaction between local moments and itinér- 
ant électrons is typically much smaller than the itinérant 
électron bandwidth. Under thèse conditions, the local 
moments tend to order into a helical state, whose optimal 
pitch K is determined by the maximum of the itinérant 
électron spin susceptibility x(q) (^he Ruderman-Kittel- 
Kasuya-Yosida mechanism,^^). As long as the supercon- 
ducting cohérence length ^ is longer than the period of 
the hélix, ^ 1, superconductivity is only weakly af- 
fected. 



II. MODEL 

Our analysis is based upon the mean-field supercon- 
ducting Kondo lattice model. 



(1) 



which has been applied with success for th eoretical de- 
scription of the bulk properties of HM^^^^. This Hamil- 
tonian describes électrons hopping between lattice sites 
and interacting with classical magenetic order pa- 
rameter S^. Here tij is the intersite hopping, J is the 
exchange interaction constant, a = (a^, a^, a^) is the 
vector of Pauli matrices, and Qq, is the operator of élec- 
tron annihilation on site i with spin a. For the helical 
state we choose S(r) = {cos Kz^sinKz^O). Even in the 
présence of helical magnetism the paring ampl itude A in 
the singlet channel remains essentially uniformP^^. 

The Hamiltonian ^ is position dépendent; however, 
the explicit z-dependence can be eliminated by perform- 
ing a gauge transformation on électrons, Q^e*^^/^ 
and Cil ~^ Qi^"*^^/^, which leaves the superconducting 
term invariant. This transformation changes the kinetic 
energy, which in momentum space becomes 

Here ë{p) = (ep-K/2+^p+K/2)/2 and h^{p) = (£p+K/2- 
£p_K/2)/2. If we assume here that the crystal structure 
is centrosymmetric; then £ is symmetric with respect to 
ail components of p and K, while h is antisymmetric with 



respect to Pz and and symmetric with respect to Px^y. 
The complète Hamiltonian is 

P 

-^Cpt^-Pi - ^*C-piCpt- (2) 
It describes superconducting électrons in a momentum- 
dependent exchange field, h = (J, O,/^^)- More specif- 
ically, the normal part is équivalent to an equal mix- 
ture of the Rashba and the Dresselhaus spin-orbit inter- 
actions, while J plays a rôle of an external Zeeman field. 
For J = there are two Fermi surfaces intersecting at 
Pz = ^ (where hz = 0). Finite J splits this degener- 
acy. The local spin quantization axis is rotated by angle 
0(p) = id.-n~^ [J /h zip)] relative to the z axis in the x-z 
plane. 

It is useful to transform the Hamiltonian to the ba- 
sis that is locally aligned with the momentum-dependent 
field h, 



Cpt 

Cpi 



(p)/2 



^Pt 

dpi 



(3) 



In this basis the Hamiltonian becomes. 



H = ^[è(p) - |h(p)|cT^^ 



1 . ,^(p) - (/.(-p) 



sm 



P/3 



(p) 



The first two lines of Eq. Q represent two decoupled su- 
perconducting bands with the spin enslaved to the local 
h(p) direction. For |A| <C max(J, the superconduct- 
ing interband coupling [the third line of Eq. Q] provides 
only a small perturbation, which can be neglected. Then, 
within a given band, the superconductivity corresponds 
to an equal spin pairing with an effective p-wave order 
parameter amplitude /i2(p)A/[2|h(p)|]. The transforma- 
tion from s-wave in the original model, to p-wave upon 
projection, is an inévitable conséquence of superconduc- 
tivity needing to conform to the lifted spin degeneracy. 
The same transformation can be applied to other com- 
binations of magnetic fields, exchange interactions, or 
spin-orbital effects. The effective order parameter will 
dépend on the détails: for instance, for 2D Rashba su- 
perconductors, the induced effective superconductivity is 
a chiral p^^ ± ip^ due to the diffé rent h(p) configura- 
tion around the Fermi surface^^^^. It is important that 
the effective exchange field acting on électrons, h, be 
momentum-dependent. Otherwise, s-wave superconduc- 
tivity projected onto either spin-up or spin-down Fermi 
surface would give zéro, as is the case for a uniform fer- 
romagnet, h = (J, 0,0)P2l. For HMS hz{p) vanishes when 
Pz = ^ and a line of nodes in th e superconducting gap 
appears at the Fermi surfac^ï^ll^. 
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III. ORIGIN OF SURFACE ZERO MODES 



J=0.25, K=1 



J=0.5, K=0.5 



In the idealized ID case (or for any fixed pair and 
Py) the Hamiltonian ([2| is unitari ly équival ent to the 
model of ID Rashba superconductoj^- ^ ^^ ^^^^^^l Since in a 
range of parameters this model gives isolated Majorana 
edge modes, from this mapping it immediately follows 
that for certain parameter range there will be Majorana 
modes on the surface of HMS. However, a more gênerai 
argument for the existence of the Majorana modes fol- 
lows from a direct solution of the Bogoliubov-de Gennes 
(BdG) équations for spinless fermions in the experimen- 
tally important régime of weak superconductivity. For 
convenience, we keep the surface perpendicular to the 
z axis, but take an arbitrary non-chiral p-wave order 
parameter, A(p). If the line of fixed {px^Py) crosses 
the Fermi surface in two points, pi^ = {Px^Py^PzR) and 
PL = {PxiPyiPzL)^ the quasiclassical BdG équations in 
the vicinity of thèse points read^^ 



- iVRÔzUR + /^rVr = Eur 
i^RdzVR + ArUr = Evr 

for "right-movers" with velocity ur = de /dp 



(5) 
(6) 

and 

similarly for the "left-movers" (the detailed dérivation is 
in the Appendix [Â]) . If the sign of vrAr is the same as 
of ulAl^ e.g., positive, then one can construct a zéro 
energy solution normalizable for z > that vanishes at 
z = as 



(7) 

Hermitian conjugation relates the zéro modes at (px^Py) 
and {—Px^—Py)' Therefore, their linear superpositions 
are the canonical Majorana zéro modes (see Appendix [b| 
for détails). If, on the other hand, urArUlAl < 0, no 
normalizable solution that satisfies boundary conditions 
can be constructed at ^ = and thus there are no zéro 
energy Majorana modes. In addition, if there is an even 
number of pairs of points where the line of fixed (px^Py) 
crosses the Fermi surface, then, generically, instead of 
Majorana modes there will be complex fermions at finite 
énergies. 

The situation is illustrated in Fig. [î] where we used 
the band structure of a helical magnet on a cubic lat- 
tice with the nearest neighbor hopping of unit strength 
and lattice constant taken as a unit of length [^(p) = 
— 2(cospaj + cos Py + cos Pz)]. lu paucls (a) and (b) we 
show a séries of energy-dependent Py = cuts of the 
Fermi surfaces for the cases J < K/2 and J > re- 
spectively (the Fermi surface topologies at low énergies 
differ in thèse cases). In Fig. [î]^, the dashed lines in- 
dicate the lines of constant Px and Py, which can cross 
either one or bot h Fermi surfaces. The case shown by the 
upper dashed line satisfies the conditions for existence of 
undoubled Majoranas: there are only two crossing points 
and between them, both velocities in the z direction, and 
the Pz-wdive order parameter, change sign. On the other 
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FIG. 1: (Color online) 2D cuts at = through the Fermi 
surfaces of a helical magnet at différent énergies: (a) for 
the case J < K/2 and (b) for J > K/2. The outer/inner 
(blue/red) line corresponds to the lower/upper energy band. 
(c) The arrows indicate the spin polarizations of momentum 
eigenstates on the Fermi surfaces. The dashed lines parallel 
to the pz axis can have either two (upper line) or four cross- 
ings (lower line) with the Fermi surfaces. The former leads to 
undoubled Majorana surface modes, while the latter does not. 
Panel (d) illustrâtes the case of arbitrary orientation of the 
helical axis, and thus also of the p-wave superconductivity, 
relative to the surface normal (z). The undoubled Majorana 
modes exist in the shaded interval of momenta parallel to the 
surface: there, both the order parameter and the z-component 
of Fermi velocity change sign upon reflection from the inter- 
face. 



hand, for the lower dashed line, there are crossings on 
both Fermi surfaces. In this case, any finite interband 
coupling will lead to fusion of the two Majorana modes 
belonging to the individual bands into a finite energy 
fermion. In Fig. [î]i we illustrate a more gênerai case of 
alignment of the axis of the non-chiral p-wave supercon- 
ducting order parameter relative to the surface. In this 
case, the Majorana modes will exist within the shaded 
range of momenta Px^Py Interestingly, such finite con- 
centration of Majorana zéro modes has also been pre- 
dicted to exist in some time-re versai- invariant systemJ^ 
the results presented here indicate that the présence or 
absence of time-reversal invariance is not the determin- 
ing factor for the présence or absence of a large density 
of zero-energy surface states. 
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IV. NUMERICAL RESULTS 

To verify the validity of the above considérations, we 
have performed numerical simulations of the HMS Eq. 
^ in two dimensions on a square lattice. We consider 
the System with a fixed number of sites in the z direc- 
tion and with translational invariance in the x direction. 
The hehcal direction K is chosen either along z oi x 
(the results for other orientations of K are presented in 
Appendix [C|). From the above argument based on the 
change of sign of velocity and order parameter upon re- 
flection, we expect to obtain Majorana modes in the K| |z 
but not in the K ± i case. It is convenient to define 
= cos~^(cosiï'/2 + /i ± J/2). For weak supercon- 
ductivity, |A| << |J|, when both bands are populated, 
p± are real and the Majorana fermions exist in the in- 
terval ( — |p+ 1 , — \p- |p+ 1 , \p- \ ) . When only one band 
is populated and p+ is real (the Fermi surface is simply 
connected) then the Majorana modes exist for px in the 
range (-p+,p_). 

In Fig. [2]we présent the results for the electronic spec- 
tra for a set of parameters representing large and small 
A/J and différent chemical potentials. The spectra are 
doubled as they include the Bogoliubov redundancy, and 
only the non-negative energy states should be considered. 
We see that in an interval of Px that approximately coïn- 
cides with the interval (p_,p+) indicated by the vertical 
yellow lines, there remains only one low energy mode. We 
verified that the energy of the mode vanishes exponen- 
tially quickly with the increasing distance between the 
edges and that the probability amplitude is equally split 
between the two edges of the sample, as expected for 
Majorana states. As can be seen from Fig. [2^, p± give 
a good estimate of the momentum range of Majorana 
modes even for rather large A/J = 0.4, and the agree- 
ment becomes excellent for A/J = 0.2 (Fig. When 
the helical direction is aligned with the edge, K||x, there 
are no zero-energy Majorana modes. The low energy 
states near p^ = are due to the gapless nature of the 
effective p^^-wave order parameter in this case. Numerical 
results for intermediate angles between the helical direc- 
tion and the sample surface are presented in Appendix [C) 
The flat section with zéro modes appears in a finite range 
of angles, consistent with the gênerai argument based on 
the BdG considération of the previous section. Numeri- 
cally we can also address the régime of strong pairing. As 
expected, when pairing can overwhelm the magnetic ex- 
change gap, the System undergoes a quantum phase tran- 
sition into a fully gapped state adiabatically connected 
to the trivial s-wave superconductor. In this régime, nat- 
urally, there are no longer any Majorana surface states, 
Fig. ^. 

V. EFFECTS OF CONFINEMENT 

For quantum information applications it is important 
to be able to isolate a small number of Majorana zéro 



modes. This can be done by shaping sample into a wire. 
Suppose that a 2D wire is confined in the x direction 
when the helical axis is pointing along z. When confined 
to a cylinder of perimeter L, the allowed values of the 
momentum along the edge are Px — 2iïn/L with n any 
integer. Since the modes with n corne in pairs and 
n = is unique, the odd number of Majorana modes at 
a given edge persists as long as there is only one (non- 
superconducting) Fermi surface. As soon as the the up- 
per energy band starts to be populated, the number of 
Majorana modes becomes even. In the case of a ribbon 
(the width in the x direction is much smaller than in 
the helical - z - direction), one can first construct trans- 
verse modes ilJn^{x/Wx)e'^P^^ from the superpositions of 
the states with Px = iTm^^/Wa^, where Ux are positive 
integers. Thèse states can be easily made to satisfy the 
zéro boundary conditions along the x axis since they have 
the same spin quantization axis. Since différent trans- 
verse modes correspond now to différent n^; > 0, even 
when both high and low energy bands cross the Fermi 
surface the number of Majorana surface modes per edge, 
Nm ^ (p+ — P-)W/7T, can be either even or odd. In 
fact, with increasing chemical potential, Nm decreases 
since the distance between the Fermi surfaces of the two 
bands shrinks. Assuming parabolic électron dispersion, 
Nm ^ J / {pF^Px)i i-e. Nm is given by the ratio of the 
exchange coupling to the ID single particle level spacing. 
Consequently, even for large total number of transverse 
modes ~ Pf/^Px^ Nm can be ~ 1. For comparison, in 
the case of a 3D bar with finite dimensions in both x and 
y directions, Nm ^ 2mJ/{ôpxSpy). Even in this case, if 
Wx ^ Wy, Nm will be either odd or even as a function 
of System parameters {/j^^Wx^y^ J). The dependence of 
the parity of the Majorana mode number on the System 
parameters, such as the transverse dimensions, resembles 
the situation occurring in the case of Rashba supercon- 
ductor s, ^^lought there the dependence is somewhat more 
intricatd^^. 



VI. SUMMARY 

We have shown that idéal surfaces of HMS host a finite 
density of zéro energy Majorana modes. When confined 
to wires, the number of Majorana modes can be reduced 
to a small number. The relatively large relevant energy 
scales - with critical températures for superconductivity 
reaching up to 10 K and exchange interaction of about 
100 K, HMS offer attractive conditions for the investiga- 
tions of the Majorana physics. 
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FIG. 2: (Color online) Results of numerical simulations of 
model Eq. ([î]) on a square lattice with nearest neighbor hop- 
ping. The System is periodic in the x direction and has open 
boundaries in the z direction. The size of the System in the z 
direction is 150 sites. The helical wavevector K is either along 
z (left column) or along x (right column). The horizontal axes 
are the momentum along the edge (the lattice constant is the 
unit of length); the vertical axes are the energy in units of 
intersite hopping. Case (a) corresponds to the intermediate 
strength superconductivity, A/J = 0.4 (J, K parameters from 
Fig. [ÎJl, the third pair of Fermi surfaces from the top). The 
range of momenta > where unpaired Majorana modes 
are expected from the small-A considération is shown by ver- 
tical yellow lines. For K\\x there are no Majorana modes, 
as expected, since the sign of the order parameter does not 
change under reflection from the interface. The near-zero en- 
ergy modes at p^c = are due to the gapless superconductiv- 
ity in the bulk. In case (b) the amplitude of superconducting 
pairing has overcome the band splitting induced by magnetic 
ordering. A fuU superconducting gap opens and there are no 
zéro energy modes for either orientation of K relative to the 
sample boundary. Case (c) corresponds to weak supercon- 
ductivity, A/J = 0.2, for the parameters in Fig. (the top 
pair of Fermi surfaces). The Majorana modes exist in the 
expected range of momenta px - marked by vertical (yellow) 
lines - as discussed in the main text. 



auspices of the National Nuclear Security Administra- 
tion of the U.S. Department of Energy at Los Alamos 
National Laboratory under Contract No. DE-AC52- 
06NA25396 and supported by the LANL/LDRD Pro- 
gram. AFM gratefully acknowledge financial supports 
from the NCCR MaNEP, QSIT and from the Swiss Na- 
tional Science Foundation. 



Appendix A: Dérivation of the Bogoliubov-de 
Gennes (BdG) équations for spinless fermions 

In this section, for the sake of completeness, we dérive 
the BdG équations (5) and (6) of the main text. For sim- 
plicity of présentation, we will assume that the Fermions 
have parabolic dispersion, and that the superconducting 
order parameter is a non-chiral p-wave with the lobes 
oriented along a unit vector h. Then the Hamiltonian in 
the real space représentation is 



H 



dr 



(Al) 



^iailj\d-n)ilj'^ ^ia*ilj{d-n)ilj , (A2) 

where ip and ip^ are the fermion filed annihilation and 
création operators, {ip^ {r)^ijj{r^)} = J(r — r'), m is the 
fermion mass, d = {dx^dy^dz)^ and a détermines the 
amplitude of the order parameter, which can smoothly 
vary in space. The Bogoliubov quasiparticles 7 are the 
eigenmodes of this Hamiltonian, i.e.. 



(A3) 



They can be expressed in terms of the fermion field op- 
erators as 



(A4) 



with u(y) and '^(r) being scalar functions. Performing 
commutation according to Eq. ( A3 ) and collecting coeffi- 
cients of ip and \p\ we obtain the following équations for 



2m 
2m 



In the vicinity of a particular Fermi point p^, we can 
factor the fast oscillating and the slow parts. 



fiju- 2ia*{d • n)v = Eu, (A5) 
- /il] V — 2ia{d • n)u = Ev. (A6) 



(A7) 



and for the slow part, finally, the BdG équations are 

— ivF ■ du + 2a* (pf • n)v = Eu, (A8) 
iiyp • dv + 2a{pF • n)u = Ev. (A9) 

The Fermi velocity, for arbitrary non-interacting disper- 
sion Ep is i/F — d£p/dY>' For non-chiral order parameter, 
and for the situation when the translational invariance is 
broken in the z direction only (e.g. by the présence of a 
surface), thèse équations are precisely Eqs. (5) and (6) of 
the main text. 
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Appendix B: Zéro energy solutions of the BdG 
équations in the présence of a surface 

We are interested in the existence of the zéro energy 
quasiparticle states locahzed near the z = surface. We 
wih attempt to construct thèse solutions as a combina- 
tion of incoming and outgoing plane waves that corre- 
spond to the same values of Px and Py^ in the vicinity 
of either pR = {Px.Py.PzR) or = {Px.Pv^Pzl) Fermi 
points. Near a particular Fermi point, we are looking for 
the solution of 

- iud^u + = 0, (Bl) 
iud^v + Aii = 0. (B2) 

(Here we dropped the irrelevant global phase of the order 
parameter A and tilde over (u^v)). Thèse équations can 
be easily integrated to find that 

u±iv = A^e^^'/r (B3) 



This is an operator that créâtes a zéro energy quasipar- 
ticle near an infinité surface for a given value of {px^Py)- 
Despite being a non-degenrate zero-energy solution, this 
state does not completely satisfy the criteria for Majo- 
rana zéro mode, i.e., 7 7^ 7^. However, for a p-wave 
superconductor with a electronic dispersion that is sym- 
metric with respect to p — p, for every solution 7 
at a given (px^Py), there is a corresponding solution 7' 
at {—Px^—Py)' They form a complex conjugate pair, 
7''' = 7^ The "canonical" Majorana zéro modes can be 
constructed out of them as 7 + 7' and 2(7 — 7'). For an 
idéal surface, there is no mixing between différent (px^Py) 
states and consequently for each allowed {px^Py) point 
there is only half of a complex fermion mode. 
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The solution has to be normalizable ai z > 0. There- 
fore, for a given sign of A/z^, e.g. positive, = 0, and 
consequently u = —iv. Including now also the fast os- 
cillating part, the individual solutions near each of the 
Fermi points are proportional to 



Note that as long as the sign of A/up is the same for both 
Fermi points, the spinor part of the solutions is the same. 
We can therefore take a linear combinations of them such 
that the full wave function vanishes at the boudnary, z = 
0. The corresponding Bogoliubov quasiparticle is 



(B5) 

I 

Appendix C: Numerical resuit s for a gênerai 
orientation of the axis of p-wave order parameter 
relative to the sample surface 



In Fig. |3] we demonstrate that the zéro energy Majo- 
rana modes exist for a wide range of the relative orienta- 
tions of the surface normal and the order parameter axis, 
as has been argued in the main text based on the analy- 
sis of the BdG équations. This proves that the Majorana 
modes generically appear on the surfaces of non-chiral 
p-wave superconductors, as long as the translational in- 
variance in the plane of the surface is preserved. 
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FIG. 3: (Color online) Results of numerical simulations of 
model Eq. (1) of the main text on a square lattice with nearest 
neighbor hopping. The parameters are the same as in Fig. 
2a, J 0.25, K = 1, M = -3.5, A 0.1. The System is 
periodic in the x direction and has open boundaries in the z 
direction. The horizontal axes are the moment um along the 
edge (the lattice constant is the unit of length); the vertical 
axes are the energy in units of intersite hopping. The size 
of the System in the z direction is 150 sites. Différent plots 
correspond to varions angles between the helical wavevector 
K and the i-axis [K is assumed to be lying in the x-z plane). 
The horizontal axes are the momentum along the edge (the 
lattice constant is the unit of length); the vertical axes are 
the energy in units of intersite hopping. Plots for ^ = 
and — ^/2 are the same as in Fig. 2a. For intermediate 
angles we see that the range of momenta where zero-energy 
Majorana modes are présent, changes gradually. 
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